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TECHNICAL NOTE 3263 


LIFT AND MOMENT EQUATIONS FOR OSCILLATING ATRFOILS 


By Alexander Mendelson and Robert W. Carroll 


SUMMARY 


Exact equations are derived for the oscillatory aerodynamic forces 
acting in an unstaggered cascade of airfoils fluttering in potential 
flow. Aerodynamic coefficients similar to those of the isolated airfoil 
are obtained as functions of the cascade geometry and the phasing between 
successive blades; the phasings considered are zero, 90°, and 180°. It 
is shown that 90° is a special case of 180° phasing. These aerodynamic 
coefficients are plotted for the special case when all the airfoils are 
vibrating in bending in phase (360° phasing). It is shown that the ef- 
fect of cascading for this case is to reduce greatly the aerodynamic 


damping. 


INTRODUCTION 


The flutter of airfoils in a cascade has until recently been pri- 
marily of academic interest. However, the widespread use of compressors 
and turbines in current aircraft power plants has given the problem sig- 
nificance. Compressor blades, in particular, are susceptible to vibra- 
tions, and some of these vibrations have been attributed to flutter. 


The problem of the flutter of a compressor or turbine blade differs 
from that of an isolated airplane wing in at least two ways. It is ne- 
cessary to consider, first, the effect of centrifugal force; and, second, 
the effect of cascading. The effect of centrifugal force can be taken 
into account with sufficient accuracy by applying the appropriate centri- 
fugal force correction to the fundamental bending frequency (ref. 1). 

The effect of cascading is much more difficult to evaluate, however, 
since it is necessary to take into account the interference effect be- 
tween the blades, which obviously depends on the cascade geometry. Two 
new geometric variables must therefore be introduced; namely, the spacing 
between bladeg and the stagger angle. In addition, since flutter is a 
time-varying phenomenon, another parameter must be introduced to take 
account of the phasing between the motions of the different blades of 
the cascade. | 
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‘The problem of flutter of compressor blades is first simplified by 
assuming an infinite cascade of airfoils. А first step in a flutter 
analysis of such a cascade of airfoils is to determine the oscillatory 
aerodynamic forces and moments acting on the cascade. Most recorded 
cases of compressor blade flutter indicate the occurrence of such flut- 
ter at high aerodynamic loading, where the blade stalls and flow separa- 
tion occurs, and it would be desirable to solve the problem for this 
case. However, no general methods for calculating aerodynamic forces in 
nonpotential flow are available. It is therefore necessary, as a first 
approach, to consider the case in potential flow at low angles of attack. 
The effect of flow separation at stall can then be taken into account 
separately, for instance, by the introduction of aerodynamic time lags 
(refs. 2 and.3) or other mechanisms which may prove useful. The object 
of this paper is to present solutions for several special cases of the 
oscillatory aerodynamic forces and moments acting on an infinite cascade 
of airfoils in potential flow. ` 


The first derivation of oscillatory forces in a cascade was made in 
reference 4. The effect of the wind-tumel walls on a fluttering isolated 
airfoil was determined approximately. This is equivalent to the special 
case of an infinite cascade without stagger, with adjacent blades being 
180° out of phase. The integral equation for the problem was set up and 
an approximate solution obtained by replacing the kernel with a simple 
polynomial. The results are not applicable for spacing-to-chord ratios 
of less than 1. The same problem was solved rigorously in reference 5 
and in reference 6 by different methods. The results are obtained in the 
form of doubly infinite series of Jacobian elliptic functions. 


Another special case was treated in reference 7. A cascade with 
stagger with all the blades vibrating in phase is considered. The method 
is similar to that used in reference 8 for the isolated airfoil. The 
form in which the final results are presented cannot be easily used for 
numerical calculations. More recently, the general integral equation for 
a cascade with stagger and prescribed phasing was set up by Sisto (ref. 
9). Approximate numerical solutions were then obtained for several cases 
of zero stagger. The method used is similar to that of reference 4. 


The present paper attempts to fill some of the gaps left by the pre- 
viously mentioned investigations. A solution is obtained for the case of 
zero stagger with the blades either 180° out of phase (considered in ref. 
6), 90° out of phase, or in phase. The last two cases can be handled 
approximately by the method of reference 9. The present paper presents 
an exact solution. Furthermore, the results are obtained in a form which 
allows а solution in a simple straightforward manner on a désk calculator, 
the final results being presented in the form of aerodynamic coefficients 
which can be calculated in a stepwise manner by the use of a set of recur- 
rence formulas. 
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AERODYNAMIC LIFT AND MOMENT 


The analysis of the oscillatory forces and moments acting on a cas- 
cade of airfoils is made using the classical methods. In the beginning, 
the cascade is assumed to be of infinite extent with arbitrary spacing 
and stagger, as shown in figure 1. The problem is later specialized to 
the case of zero stagger. The airfoils are assumed to be thin and per- 
forming small oscillations in a potential, incompressible, ideal air- 
stream. The airfoils and their wakes are replaced by surfaces of discon- 
tinuity (vortex sheets), the interaction between the vortices being 
neglected. Each airfoil is performing both bending and torsional oscil- : 
lations, two adjacent airfoils being out of phase by a prescribed amount. 


Under these conditions Euler's equations of motion are first linear- 
ized and then solved for the pressure distribution over the airfoils. 
The Biot-Savart theorem giving the induced velocity at any point on the 
airfoil due to the vortex field must also be used to obtain a solution. 
Once the pressure distribution is known the aerodynamic lift and moment 
acting on the airfoil can be obtained by integration. The complete solu- 
tion is given in appendix B. The final equations and results will be 
given here. 


General case. - The aerodynamic lift and moment acting on an unstag- 


gered cascade as obtained in the manner outlined can be expressed in a 
form similar to that given in reference 10 for the isolated airfoil. 


L = προς In Et ls - (5. a 3 


M = pb o (|M, -(& + ΠΣ Ë = (2+ Je, +) + (F + a) In α 


where (all symbols are defined in appendix A) 


а, elastic ах1в position as fraction of semichord Ъ, measured from 
midchord, positive toward trailing edge 


b semichord 
h bending deflection, positive downward 


α torsional deflection about elastic axis, positive for increasing 
angle of attack 


p mass density of air 


(Ὁ circular flutter frequency 
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The aerodynamic coefficients Tes x M M, are functions of the 
following parameters: 


k reduced frequency, wb/V 
21m phase angle between adjacent blades; m = О, 1/2, 1/4 
À cascade geometric parameter given by À = 1/8 


For the isolated airfoil these coefficients are functions of only the 
reduced frequency К. For convenience, the following additional param- 
eters are used: 


_ À 
g = coth 5 
ε = coth À 
x = 641 

“© E-l 
т = 2-1 

τ 

= ik 

x = co ο 


where 
8 spacing in units of semichord 
V free-stream velocity 


The following scheme, as outlined in appendix C, is now used to calculate 
the aerodynamic coefficients: Let 


κα... 
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ΡΝ 2: cj+l+21 * ©3-1-21. А T 
E .Jtivei '^^d-l-e6i | > 
J ΞΟ — (zui)! -J 
es 
атара + 8321-21 
е, 2i = e J 20 
i=0 (2i+1)o “J 


f, = с ... 
j J- 17 111. 


С, = (2) (ез-1 Р AC 4 + C 443)£3 


j 
Е, = (-1) (611 + 2e, + e 443), 
Ё, = (-1)3 (à, 1 T ed. + δα) χι dy 2t 3, 1-21! 
б, = (2) (а, + 24, + а) b° (ej οι + 941-21! 
і=0 


Ay = 
2341 
Азы. = 234 А; 


r(z) = gamma function of argument z 


(x,j) = х(х+1)(х+2) + - + (x+j-l) = I 


» l/2, 1)ix, і 
F(1/2, x; x + 1/2 + q; 1/τ}Ξ > x 41/2 + q,i)(1,1 
10. 


is the ordinary hypergeometric function. 
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The aerodynamic coefficients can now be obtained as follows: 


| P 
1 1 ай 
i o tae S 


上 


_ fl mW ` j 
πο, 


κ 


3 


= +W - SL Y+ =L - P. + 4р, +P +2 (ος p) -È pD 
М = 2 αν ως us 


The geometric parameters À, €, 0,7, T are functions of the spac- 
ing and have been tabulated in table I for different spacings. 


The quantities a; and b; are functions of only the phasing m. 
The first 20 values of bj bave been tabulated in table IL for values 
of m equal to O and 1/2. It is to be noted that for m - O, b; is 
equal to 1 for all J. 


In the equations for cj， dj» and ejs the index J never takes on 
negative values. The quantities ET d. 1 e. must always be obtained 


by calculating C4, dy, and es. 


The Ру, Po, Ps, Py, and Ps are functions of the cascade geometric 


parameter А and the phasing between blades m. The functions С and 

D are functions of A, m, and the reduced frequency К. It is shown in 
appendix B that the case m = 1/4 (90° phasing) reduces to the case 

m = 1/2 (180° phasing) with twice the spacing. For the case of an iso- 

lated airfoil, the functions reduce as follows: 


в» 299 P; + Ó | 
À +0 P, + 1/2 
Ру» 2 Ps + 1/8 
Р» 1 
C(A,m,k) + C(k) Theodorsen's function 


2 
k 


оак) > - cf) - Ж |: - с(к)] + À 
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The coefficients Τη, L,, M, and M, then reduce identically to those 
of the isolated airfoil as given in reference 10. 


The values of À; are independent of geometry or phasing. The 
first 30 values are tabulated in table IIT. 


Special case, m = O: For the case of all the blades vibrating in 
phase (2-0), & few simplifications can be made. 


a 
g -J 
ο. = σ 
J g^-1 
a = 
f, = 2 κα 
(c-1) 
Ῥ = 2 
1 7 Xe-1) 
£ ef 


-1/2 
2 1 {5 (2 x F (1/2, x; x+1/2; 1/+) 
€”? — 1 x+ 1/2 
C = K Jn ——— + — ο —— 
2 1/2 
ε T F(-1/2, x; x + 1/2; 1/*) 
For this case, therefore, the bending coefficient IL, is obtained 

in closed form. The preceding results -can be obtained by summing the 


series for the appropriate quantities or by integrating 14, Ie, Ig, and 
17 of appendix B directly for m = 0. 


APPLICATION AND RESULTS 


As an example, the case in which all the blades are vibrating in 
phase (m=O) will be considered. Experimental data on compressors indi- 
cate that the only important type of vibration occurring is one of pure 
bending. For this case, then, the only aerodynamic coefficient that 
need be considered is Ly: 
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The functions Ру and Р. are plotted in figure 2 against the 


cascade geometry parameter À. The real and imaginary parts of the cir- 
culation function C are plotted in figure 3. It is to be noted that 

C is a function of both the spacing and the reduced frequency, whereas 
the P functions are independent of the reduced frequency. 


in figure 3, over the range of reduced frequencies from zero to 0.6 
and spacings s = 1.67 to 3.33 (this corresponds to solidities of 0.83 
to 1.67), the real part of the circulation function C is independent 
of reduced frequency for practical purposes, whereas the imaginary part 
of C varies linearly with reduced frequency. 


With the P functions and the circulation function C available, 
the aerodynamic bending coefficient Lp, was computed and plotted in 


figure 4. The values for the isolated airfoil are also plotted on this 
figure. It is to be noted that L, (as well as the other coefficients) 
is complex. The real part of Ig 15 in phase with the displacement and 
can do no work; the imaginary part is in phase with the velocity and cor- 
responds to the damping component and can do work. Whether this system 
is stable therefore depends on the sign of the out-of-phase component or 
imaginary part of Ly The real part of IL, is practically independent 
of reduced frequency and varies only slightly with spacing in the range 
considered. The imaginary part of Ту, however, varies appreciably 
both with reduced frequency and spacing. This is further illustrated 

in figure 5, where the coefficient L, has been plotted against spac- 
ing for a reduced frequency of 0.4. The results of figure 4 agree 

well with those of reference 9 for а spacing в of 2.0, which is the 
only one considered there. ‘Some of the data from that reference are 
plotted in the figure. 


Since the imaginary part, or out-of-phase component of L), corres- 


ponds to the aerodynamic damping, figure 4 shows clearly that the effect 
of cascading is to reduce greatly the aerodynamic damping when the blades 
are vibrating in phase. For example, at a reduced frequency Kk of 0.4 
and spacing s of 2 (solidity of 1), the aerodynamic damping is approxi- 
mately one-half the value for the isolated airfoil. At a reduced fre- 
quency of 0.1, the damping is reduced to almost one-third the isolated 
airfoil value. The higher the solidity, the lower the aerodynamic damp- 
ing becomes. However, for this particular phasing (all the blades vibrat- 
ing in phase), the damping never actually goes to zero in pure bending 
except in:the limit of infinite solidity. · 


CONCLUDING REMARKS 
The equations for the oscillatory aerodynamic forces acting in sev- 


eral cases on a cascade of airfoils in potential flow are derived. The 
aerodynamic bending coefficient is calculated and plotted for the case 
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where all the blades are vibrating in phase. For this case, the effect 
of cascading is to reduce greatly the aerodynamic damping force. 


The results of this investigation can be used in the study of stall 
flutter phenomena in a cascade in а manner similar to the use of the 
classical flutter theory of the isolated airfoil for the stall flutter 
of the isolated airfoil. Aerodynamic time lags can be introduced into 
or characteristic times (ref. 2) can be deduced from the aerodynamic lift 
and moment equations. The linear part of the blade characteristic (ref. 
11) is directly obtainable, and the nonlinear part, if known, might be 
treated by the methods of reference 2. The general usefulness of this 
approach, for both the isolated airfoil and a cascade, must, however, 
still be determined. 


Lewis Flight Propulsion Laboratory 
National Advisory Committee for Aeronautics 


Cleveland, Ohio, July 16, 1954 
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APPENDIX А 


SYMBOLS 


The following symbols are used in this report: 


λος, 
cj 34; 6) 


b 


= À,m) 

D(k, 2,m) 

"Pr t d" 
J^ 30 


< 


f(t), f(n) gfu) 
h 


i 


ЕМ. 
Το» Т6? 7? ТӘ 


elastic axis position, measured from midchord, as а 
fraction of semichord b, positive toward trailing 


edge 


coefficients in recurrence formulas given in text 


semichord 


functions defined in text 


coefficients used in evaluating integrals, given in text 


base of natural logarithms 


functions of indicated variables 


bending displacement of airfoil, positive downward 
γ-ι 


integrals defined by eqs. (B40) and (B41) 


reduced frequency, wb/V 
lift per unit span 


aerodynamic coefficients in lift equation, defined in 
text 


moment about elastic axis 


aerodynamic coefficients in moment equations, defined 
in text 


phase angle lag between any two adjacent blades, as 
fraction of 2x radians, 5/2x, appendix B 
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functions defined in text 


pressure at any point on airfoil, function of time 
amplitude of pressure difference 


pressure difference at a point on airfoil, p,' - р,', 
function of time 
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spacing between blades in units of semichord hb 
geometric parameter defined by T 

time 

local velocity component in free-stream direction 


free-stream velocity 


local velocity component perpendicular to free-stream 
direction 


coordinate in free-stream direction in units of semichord А 


coordinate perpendicular to free-stream direction іп 
units of semichord 


angular displacement of airfoil 

stagger angle | 

total circulation about airfoil 

vorticity distribution for reference airfoil 
vorticity distribution in wake of reference airfoil 
phase angle between two adjacent blades 

geometric parameter defined by z=' coth 入 


transformed variables of integration defined in appen- 
dixes C and D " 
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Subscripts: 
1L,j,n,p,q 
u 


7 


geometric parameter defined by Ξ e 
air density 


geometric parameter defined by o = coth = 


geometric parameter defined by т = 


η] 
PIS 
| 
ue 
— 


velocity potential 


circular frequency of oscillation 


summation indices 
upper surface of airfoil 


lower surface of airfoil 


The subscript n is used also to indicate the neh airfoil 


一 一 ------ ----- 一 一 一 
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APPENDIX B 


LIFT AND MOMENT EQUATIONS 
General Theory 


The analysis of the oscillating airfoil of infinite aspect ratio 
in a cascade proceeds along classical lines. The airfoil is assumed to 
be of small camber and thickness and performing infinitesimal oscilla- 
tions in an incompressible ideal fluid moving at a velocity V at 
infinity. 


The airfoil and its wake are replaced by a surface of discontinuity 
or vortex sheet of strength ү! over the airfoil and τ in the wake. 
This surface of discontinuity is assumed to lie in a horizontal plane 
parallel to the direction of flow. The vertical displacement due to the 
interaction between vortices of this sheet is meglected, this displace- 
ment being assumed small compared with the horizontal motion. The 
assumption of small perturbations to the free-stream velocity permits 
the linearization of Euler's equations and, with.the introduction of a 
velocity potential, Bernoulli's equation for nonsteady incompressible 
flow is obtained. The derivation ‘is given here in detail for completeness. 


Euler's equations for two-dimensional flow can be written 


È (viu) + (Vu) È (μα) + v 2, (йыз T: 
(B1) 


ov dv dv _ _ 1 Әр! 
oe жр 1 Әр" 


Considering only first-order terms and realizing V is a constant give 
the linearized Euler equations as 


ou ou | 1 dap! 
Sa ο 
МЕ " (B2) 
1 t 
S6 * ax 5 g 
The velocity potential ф is now introduced: 
ека οφ 
= (B3) 
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Equation (B2) can now be written 


ð (οφ ὃ {ὸφ : 
= (#)+v2 (š T- Xx 


А (B4) 
ὃν (8 Va) 
which leads directly to Bernoulli’s equation for nonsteady motion 
dP ΠΕ 
Ἔν 28+ 2р = f(t) (B5) 


The airfoil and its wake are now replaced by a vortex sheet of strength 
y! on the airfoil and yo in the wake. The difference between the 


velocities of upper and lower surfaces of the vortex sheet is therefore 
Υ' or Yy- For a point on the airfoil, therefore, 


da Әр o 
ox dx -Ί 


which upon integration gives 
X . 
mme f r' dx, (вв) 
Equations (В5) ала (В6) then give 
3 x 
αρ’ = won το εν δε [ ἵνα ax, | (B7) 


Continuity of pressure in the wake requires that Ар be zero there, and 
equation (B7) becomes, for the wake, 


: b or! 
tte f Тү ОАЕ ТОЕ (58) 


where 


1 ` 
r' = f. Υ' dx (B9) 
-l 
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The assumption of simple harmonic motion will now be made, so that 


y! = yet 
ν' = velot 
popu 
[on = Y. e t (B10) 
1ος Tost + 
Гт=ге =e J. Υ ax 
-1. 
Equation (B8) now becomes 
е | 
Y. + dk rà T. dx, = - 1КГ (B8a) 


Equation (B8a) is a simple nonhomogeneous integral equation with the 
kernel equal to 1. Its solution is : 


~ik(x-1 
ү = -iklTe G-1) 


(811) 
Equation (511) gives the vorticity distribution in the wake as a function 
of the total circulation around the airfoil. It can also be obtained 
from the condition that the total circulation around the system compris- 
ing the airfoil and the wake must equal zero. 


With the use of equation (B7) the pressure distribution over the 
airfoil can now be written in the following form: 


x 
rmx р г ax, (512) 
-1 


The Lift and moment about the elastic axis can be found by inte- 
grating the pressure distribution as follows: 
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1 l x ' 
L = -pVb J, y dx - ipVbk / ydx, dx 
“| -L  f-1 
| 5 L А | L + oX 
M = -pVb J. r(x-a)ax - ib “pVk d (x-a) T Tdx. dx 
‘Wel к ~1 -1 


Vorticity Distribution in a Cascade 


(B13) 


, Consider an infinite cascade of airfoils of chord 2, stagger angle 
B,-and spacing s as indicated in figure 1. The airfoils are replaced 
by vortex sheets, the vorticity of the nth gheet being designated by 
use of a subscript n. The Biot-Savart theorem is applied to give the 
induced velocity at a point x on the reference airfoil due to an element 
of vorticity of strength ү, located at the point (χι, ns cos β). 


=y (x - x, )dx4 


Ду = SS ————————— B14 
2x|(x-- κι) + (ns cos В) 
Measuring trom the y' axis, this can be written 
"M Y xx = ng. sin В)ах 

-dv,, = = = 1. E (B14a) 


21|(ns cos в) + (х - ху : ns sin où 


Summing over all n and integrating from -l to * change equation 
(Bl4a) to | | 


κ Е Y (x 2 x, - ns sin pax, 
v(x) = - = [us cos B + Guns sàn p)] [2 сов p саа asa p) | (B15) 
---- | ns cos B + i(x-x,-ns sin B)| ns сов 8 - i(x-x,-ns sin В) 
zl . ы ^7 5 


The assumption is now made that any two adjacent blades are out of phase 
by an angle | 


Ὁ = 25m O<m<1 


---.π----- -- кч -- w ., с... 
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This angle is constant through the cascade. The vorticity distribution 
of the n”? blade is therefore out of phase with that of the reference | 
blade at the origin by an angle 


‚б„ = Zum 
or 


тете = (B16) 
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where y is the vorticity distribution for the reference blade. Sub- 
stituting T into equation (B15) results in 


"ατα 5 us m ea ὅτ (817) 
n=- n 一 | Ir-en ~ τῷ 
-1 _ se -L se 
The infinite sums can be evaluated (see ref. 12). 
= 2mm - -A(1.-2m) (x-x4 ) Ë 
———— ον ee (B18) 
- е1В(х-х,) Lu sinh A(x-x.) 
is 
= X ip 
À = - © 
Substituting now equation {B18) into equation (B17) gives 
E -A(2m-1) (x-x, ) _ б ÀA(2m-1)(x4-x) 
mer) — mail í dx 
у =з Se 
4x ginh Ax, - x) L sinh A(x4-x) 1 
-l -1 | 
(B19) 
= WE ou 
мее 


It should be noted that although equation (B18) is not valid for 
m = O, equation (B19) is valid for all m. Equation (B19) gives the 
induced velocity at a point on the airfoil due to the complete vortex 
field of the cascade. It reduces to the well-known equation for the 
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isolated airfoil when the spacing s goes to infinity. If use is made 
of the condition that the normal induced velocity must equal the normal 
component of the velocity of the airfoil, then equation (B19) gives the 
relation between the normal component of motion of the airfoil and the 
vorticity distribution. The problem then becomes one of solving equation 
(B19) for the vorticity distribution as a function of the airfoil motion. 
Once the vorticity distribution is known, the pressure distribution can 
be calculated from equation (B12). 


The solution of the general equation (B19) will not be attempted 
here. Instead, several ii ang cases of zero stagger will be considered. 


For this case, В = 0, А = = xs. , and equation (B19) reduces to 
"Y )cosh A(1-2m)(x4-x) | 
v(x) = А 一 一 一 一 一 一 一 一- ах, (B20) 
sinh sinh cx) -x) 


Equation (B20) will now be solved for the following three values of m: 


m= 0 all the blades are in phase 
m = 5 adjacent blades are 180° out of phase 
m = = adjacent blades are 90° out of phase 


Τα 
o 
H 
B 
lI 
O 
` 


x P cosh A(x,-x) 
v(x) = Эл { y(x) šini M(x, -х) dx, 
À ( À 097x) : 
κ TZL Sinh A(x, -x) ш! 
1 


(B21) 


The second equation follows from Kelvin's theorem, since L y(x} )аху=0О. 
For m = 1/2, | -1 


= ὃς — dx, . (B22) 


Sas". -- =e _` «= = — eee mn баила отт щл eee n et -————— аа - 


— — - — — € ---ν-- MM — s m — GI cda  — + —— -Ἡ--- —s — —— — — — — ------- = 


--- ------.---- ------ mm 一 
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cosh AG -x) 
vex) ~ 2л ) sinh sinh λ(πη-κ) = 
| rey) 

К sinh 4(x,-x) 


For m= 1/4, 


dx, (B23) 
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^ The case for № = 1/4 therefore reduces to the case of m = 1/2 with 
twice the spacing. The same result can be obtained from purely physical 
considerations. Thus only equations (B21) and (B22) nave to be solved. 


In order to solve equations (B21) and (B22), a transformation of 
variables similar to that used, in reference 6 will be made. Let 


- | tanh Ax = u/e 
tanh Ax, = n/e | (B24) 
‘tanh À = 1/6 | i 


Equations (B21) and (522) can ben be written in terms of р and η, 
making use of equation (B20) , to- become 


| | RI ^ 
v(u) (ea ET [ie тар -ᾱ γ(η)(ε-η) an 
: | 


pr ie πο. 
xar K = P (e2-n2Xe-n) В" 


+ m (B25) 


Чеге m takes on the value of О ог 1/ 2. Or, more briefly, 
i | Т cn Ε(η) 
g(u) = 2x J an (B26) 
-1 


where g(u) and f(n) correspond to the appropriate parts of equation (B25). 


Equation (B26) has an explicit solution for f(n), first obtained 
in reference 13 and proven rigorously for the real domain in references 
14 and 15. With the condition f(1) finite, the solution is 
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1 | 
ο JE | ε(μ) лу 28 Se (ват) 


The solution of equation (В25) can now be written directly. 


1 
rp. Ton (em) vu) 4 Bt Gant ay 
= (exh jp l-u (еар уб η-μ 


1κ1 (ε- (en) dn; 
2x ^ ИА тї πα η-μ ii (e+n,)” (ern )*. EM ut bae 


The second integral can be partly evaluated by reversing the order of 


integration: 
qti +1 
1-р \т-н H-1 cm n 


Substituting equation (B29) into equation (B28) results in 


1 
ms р -1 
τη) 2 a [in (eza τίμ) AEE (eza) ар 
£^ n? х ltn (e-n)° l-u (εἰμ) η-μ 
p r= AL +1 Su. 
[ (stn, ) 1 1 
also 
Ze RN an (B31) 
1 А 


Equation (B30) can be integrated as indicated in equation (B31) to obtain 
the total circulation around the airfoil T. 


mm —— mm rm :. ο лыш — em ч M — + ——EFIAIWILP ` — — n n a e + τπτ  — — H 一 -一 = 一 
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2e [ Fer ^u) AE (em а à 
a 7 (ew ga M VP 1 


ik m-1. Natl (ε- ju dn 

l + = x T үка em Γ сле n E NUN 

ET Е (8-1) A (en, y 7 
(B32) 

Equations (B30) and (B32) give explicitly the vorticity distribu- 


tion over the airfoil. With the vorticity distribution known, the pres- 
sure distribution and hence the lift and moment can now be calculated. 


Lift and Moment on Airfoil in Cascade 


Consider the airfoil shown in figure 6. The upward displacement 
at any time t of a point x on the airfoil is 


y(x,t) = - h - b(x-a)a (B33) 
The induced velocity у(х) must satisfy the condition that the flow 


is everywhere tangential to the surface. If the induced velocity in the 
x-direction is small compared with free-stream velocity V, this condi- 


tion leads to ( s f y) 


š v(x) 


| (554) 
= - (h + Va - ba& + bax) 
or, in the η coordinate, 
ν(η) = -(h + Va - baa) - = a log 2 | (B35) 
The lift and moment are given by 
1 b 
L= b J 名 Ca =" / д }ак 
-1 -b 
(B36) 


1 


M = b“ ; ASCH Gea ex 


1 
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Combining equations (B13), (B30), (В52), (B35), and (В56) and per- 
forming considerable algebraic manipulation, the Lift and moment can be 
written in the following form: 


zt [ы - (244) a) 
κε κακό, 6s] ee [ж -(%+ «) зы (3 9) EM 


(B37) 
where 
Ls P (m) - = P. (A,m) C(k,A ,m) 
Е (*,m) 
L 三 5 Γη + P.(A,m) - = P (A m) 1 + οίλο, κ) - m" C(A,m,k) 
M, =(3+ + i - Р, (^, m) + 一 一 -一 一 ЕЕ rm) D(A,m,k) (B38) 


i= at. + M, - 2 L) + = L, ` P_(A,n) + [P4 0m) З] + 


[cum x) + рк) | = ^s P. (A m) D(k,A,m) 


- -- ————————————— —— —————————————————————————————— ane eS co gee ey ey ee mm -- --- 


24 NACA TN 3263 


LE 
P. (Kam) = - = 
у xÀ E 
2m) = - — I, 
X À 
P2(A,m) = - —— T 
5 элэ 5 Š 
GQ 
E 
(B39) 
£ 
ik 
1 + ік + — k^ e E. 
οπλ 
C(A,m,k) = ik 
1 + і £ ke Le 
1 + ik - 1 xê _ de κὖ QE. 
а απλό 3 
D(A,m,k) s - = 
& 
1 + i Ax ke. 16 


and the I's represent the following integrals: 


(ε--η ΠΕ 1- š ( ξ-μ - Tu a 
dq (ε-η)α FE Ἢ (ew) ϊ-μη-μ | 
| (en). [a (e, fu, ењ ἂμ à 
2 | κ (ean) ltr 1 1-р 798 επ. η-μ a 


(eu). 


III 


M 
UI 


= (en) Es ε-η k g Sh ἂμ. α 
+ i (ε -т)® 1+1 log 611. (eru)? ТО? 8 е 1-и an 
ы GT f {1ος sa ur [n au 
14 E ТА (ε-η) 1 Γη (тов =) m (etu)! l-u η-μ an ' 


= 
l 


_ (є+ (ex) m-l - 


(B40) 
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m " (es ay VER | SA es 11. any 
u (e-n)™ Vien (em Ves m^ 
= ΠΡΙΝ ΤΙ (ua | (e-n, *-4 faiti dn; 
е " -1)7 L se) (e+n,)* 1-1 11-1 
E ae s 

Ta = (any 4 fn Үш (то єз} emu [1 x 
a = —————— — 

(eg) Уп (etm)x Л 11-1 


(B41) 
It is to be voted in the preceding integrals that 
lel zi Jul <1 
lili [nf 22 


The integrals in equation (B40) are functions of only geometry and 
phasing between blades (À and m). The P functions are therefore 
functions of geometry and phasing only. The integrals in equation (B41) 
are functions also of the reduced frequency Е. The C and D func- 
tions are therefore functions of reduced frequency as well as of geome- 
try and phasing. The integrals in equation (B40) are evaluated by inte- 
gration around & clósed contour in the complex plane. The details are 
given in appendix C. The integrals in equation (B41) are evaluated by 
means of the hypergeometric. integral. The details are given in 
appendix D. 


In the limiting case of the isolated airfoil (А = О), the following 
limiting values are obtained for the various functions: 


170 ῬΡη(λ»π) > 2 P5(1,m) + 1 
Р5(А,ш) + 0 Pa (X,m) t 1/2 Ps(A,m) 1/8 
C(A,m,k) + C(k), Theodorsen's function 


k2 


D(A,m,k) + - C(k) - = (1 - c(k)) + "E 


To Sn — Y —- à re — mM cs t. з 


inm — wt ELM 
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Equations (B37) and (B38) then reduce to those for the isolated 
airfoil; Ip, Ig, M, and M, reduce identically to the coefficients 


for the isolated airfoil given in reference 10. 
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27 
APPENDIX С . 
EVALUATION OF INTEGRALS Τι TO Is 

The integrals I, to Ig are evaluated by integrating around а 
closed contour in the complex €-plane, where the mapping function re- 
lating the n-plane to the C-plane is given by 

| К: τος м 
η = a(t + 1/6) x (ci) 


By use of equation (Cl), equation (B40) can be written as double 
integrations around the closed contour C (fig. 7). Let 


T 2 _ „121 = à 
S; = є + А/6°-1 Co = € - Nef- uc 


tz = - e + Ме2-1 =- t, =-e-a/e21=-t 
5 t 4 1 


Then 


(CAE can ('- ТР (CHARTES Les a1)? 
(5- eR CS t5)" (ζ'- TOS 6'-5)(5'-1/5 


(£-c,)"-1(c-c, ym-1 EE а [estes (0-0) (0-6) (p+)? | ELS -ζκ)ίζ'-ζι) 
(ζ- σι) Ct E) (CS CT LS 5η ὑ [ee c')(c"- ш 


ret 

C 

Е" | Gr MCE)? (ery? (а-006-00) l (c-r) 
C 

„= [ 
© 


dot a - 
аст at 


gar зов (ESE ME) EN 
(Re f GG) CT σι WO ($-0)(0-0) 


(£-C3) "3 (6-04)? (en)? 3 (61-6) (6-52) 2 И! Се) (0-671 


(6-65) (6-64). а" 


(p+)? 
Ge $ (r-t Eg (EEA) 


(G'-E5)(6"-64) 


ζ-ζ yi Ç- m-1 2332 = „г. yB-l/ ni m-1 2 
E (£-63) (5 人) (p) MSS (ζι-σ)(ς ЧА (ς'-ζῃ) G шы ως (сє) (er) MUST LM ΠΡ 
(ζ'-ζς)ῆίς'-ς,)5 (Ε'-ε)(ζ'-1/τ) (€,-6') (6-0) 


(6-63) g-t)" ; (6006-04) 3] 
(c2) 
lt is to be noted that in integrating around the closed contour C, the 


integration path С! in the r-plane is actually being traversed twice. 
Bach integration around C mast therefore be πα ay 2. 


If the integrations in (C2) are to be performed, the orties of 


the various integrands must be considered. It will suffice to discuss 
one of them, since they are all of the same general type. 


———— ————  ——À — — —— —— m m 
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Consider the integral Τσι defined by 


一 


(с eet ы. (t'-ta) (6-64) ate 
se | ατα CCD CUE 
| (C3) 


The integrand of Iz) has branch points at SE Co = 1/64, ts = - 1/64, 
and Ü, = - 1 as well as simple poles at ζ and 1/ζ. The transfor- 
mation (Cl) maps every point in the n-plane into a pair of reciprocal 
points in the C-plane. That part of C which lies below the real axis 
is therefore the reciprocal of that part of C which lies above the 
real-axis, and every point which is outside the contour С has its 
reciprocal point inside the contour. The singularities of the integrand 
&ppear in reciprocal pairs, one of each pair lying outside the contour 
and one inside the contour. The poles at Ё and i/t lie on the con- 
tour (see fig. 7). 


Considering the branch points, branch cuts are now made between L, 
and ὃς and from 0, and to infinity as shown. The integrand is 


now unique and single-valued over the path of integration. Furthermore, 
it can readily be shown that since the integrand of Iz, is invariant 
under the mapping ξ' -- 1/ζ', (this evidently mst be so because of (C1)), 
the sum of the residues at Ü and 1/6 must vanish iđentically. Con- 
sequently, any other contour such as Со, equivalent to С but exclud- 


ing the points ζ and 1/t, may be used. 


The integrand of Is, is therefore analytic in the annular region 
enclosed by the curves Су and Co as indicated in the figure. It can 


therefore be expanded into a Laurent series about the origin valid in 
this annular region. If Ау is the coefficient of the 1/t' term in 


this series, then | 


In order to expand the integrand into a Laurent series about the 
origin valid in the indicated annulus,.equations (C2) will be written 
as follows: ΠΠ 


= = — = 2 
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- — d. —.. ранд е 
home ARS -ν----- анне mua - de аы — — A tomes L. —— mens . w mi - me чл LA лае τον 


m-l m-1 n-1 a-l 
1 + À 1 +Š 2 τ 1 Q- l (03 1)* 


| | 
= η t 
ee ὃν pae Cem 
i Cr NOI i G-u) 0-8) : 1 +217) ( + [Д (C + 1) 
=-= = - = ac! d 
ee aes ` 9 | ee ын: уау“ 
О 


1 m-l а pli 1 [4 2 1 ю-1 NE 1 [44 1 2 

1 + 1 + 1 - 1 - 1 - 1 - 1 + 1 + 1 + 1 

( 2) ( j 1 3) ( 3 (1 s 5) ( str} í --) log | zal =) FT (š ) dt! ac 
g т ; 


ἀζ! ag 


с92= KE VOVN 


І, = - bs og - | | | 
"Jd μα) 2+2) +š) б+ 2) 642) 0 -3)6-$)/6- Ὁ ᾱ- το) 

А 144 UG + үс ü - à JQ - š) à EC ü - = ^ - ὯΝ {r (t + 1)? 
[ES 1 p DEI d t (6 - 2) 1 ντ AE a SERIES 

с ex)t-i б +) + š) Q+ar) 642) Q-go en 

š ü + à ας ТШ de Q I 5) 4 112 ü = = Mc я gy 1: pete) (a + 如 l (ри + 1)° 
Is = - =s = — | 108 U (1 Е 5) E | 1og — аста; 

“J ещ 6X9 ΝΩΕ ет знн 

(cs) 

vbere 


L 1 YAE À 
а = кє — = - = = - = g + \/s -L = coth = 
S бо ὃς δ, 2 


Each integrand can now be obtained in the form of a Laurent series about the origin by expanding each 
bracketed quantity (noting that each of these series will converge at least within the annulus) and 
multiplying the resultant series together. In this way, for example, the coefficient of 1/t' can 
be obtained and the inner integration performed. This coefficient will, in general, be an infinite 
series in Ё. The process then has to be repeated to carry out the second integration. The final 
result will be а doubly or triply infinite series in G. The process is evidently extremely long and 
laborious. However, it can be greatly simplified by employing а get of recurrence.formulas which can 
be obtained by inspection of the series, whereby each term of a given infinite series is obtained in 
terms of a previously given series. The calculation scheme is ag follows: Let 
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m у Cj+lt2i + Cj-1-24 _ d 
jJ £L (2а) 7d 


o 


Jjd2i + d3-1-2i _ 
| Bh 
4—0 (2441 ος 


i=j 
Dj = (31) (Q4. + 24, + τς ке 2, a) 
i=j _ 


— 


Е, = (21) (65.3 + ге, + μι + < > D 


i=j 


Sg 


了 -二 
F; = (ay (41.1 = ad, i 8.341) 2 (à, 2: Ы 3, 1.81) 
1 


J-1 


(Gy = (07 ὥρα + 2, + ἄγη) > (63-21 + 31.21) 


130 
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APPENDIX D 


EVALUATION OF INTEGRALS lg TO Ig 


Let 
у= if 1/8 e ) 
ετη 3 
со 
Je em š 
r 1 η. 
Yi = £T (D1) 
l 
. ‚1/2 = etl 
£ — &-1 
т = v 
^ 
q = τ-ι 
(0879. | 


The order of integration is reversed and the logs are S panden in order 
that equations (B4) can be written as follows: 


l-x-m i 1 
тт ^ х-1 ME үн ° -1/2 1/2 -m 
Te eem Д (1-v. ie ν (1-v) “(1-Ty ) 
O O 
(1-т'у) À dv ду 
al 
l-%-m ϱ 上 1 
zi j-1/2,, 31/2 
ο πα ê >` q | хе ας m /2 ν (1-v) 
т ш „2 Ут ды т 
(1-Ty) ` T(i-T!y)* -l ἂν дуу + 5 = т" log T 
l-x-m l 
y e 
sss o. PE was Е i үн ' „1+5-1/2 
8 2e 
» » O 


2 2 
(a-v) 4(1-1y ) ee ἂν дуу + L; log ç - i Ig log τ (D2) 
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The inner integral is the same for all three, except for the expo- 
nent of v. This integral can be evaluated in terms of Appell’s hyper- 
geometric function in two variables of the first kind (refs. 16 and 17). 


© ple 1/2 m -1 (n 1/2)r(3/2) 
| 1 к (1-v) ° (2-ту) (1-ту) ἂν ου 


Fi (n31/25m,1;n42 ST pt) (D3) 


Г (n) Gamma function of argument n. 


The hypergeometric function F) is defined in terms of a double 


infinite sum in Т and "n. 
€ ee 
Ру (a+l/2;m,1;n+2;T, T) = ), » туул TUE P(r)? (n4) 
p=0 q=0 nte,ptgJp:. 


where 


T inc 


lin (ps) 


(n,p) = n(n+1)(n+2) © + -(п+р-1) = 
(n,0) = 1; (1,n) = n! 


F4 can also be expressed &s & single infinite gum of ordinary hyper- 
geometric functions. 


. Se | 
ERIC) T’ 
F (n+1/2;m,1;n+2;T,T') ὃς 2 nt2,p)p! T? F(-p,1514m; 1 -~ =) 


(26) 


1 
where #(-v,1514m; 1 - z) is a, poLynomial given by 
oD 


SERT ης > = ( - 工 -19 
F(-p,1jlim 1 -g-) = Img U. s) 
q=0 
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Upon substituting back into equations (D2) ; the following integral re- 
mains to be evaluated: 


Γ уау) ; үм 2-9 Ἢ 


O 
_ LT +g T (1/2 F(q Е 1/2, x+q; X+q + 1/2; 1/х) (рв) 
Г(хна+1 /2 2 ” ᾿ 
Г{х+а /Г(1/2) aa 
= келд T Ε(/2, x; x+1/2+q; 1/τ) 


(ref. 18) 


Equations (D2) can now be written as follows: 


l-x-m - 
2 М 
T c LO JP (1/2 1 dtp 
一 -一 一 一 一 -一 ) | | = À p T + 2AL 
| + ud jtp 6 
x Jel O ? P 


.H 
II 
à 
H> 
£n 
g 


l-X-m ә oe em 
2 4 
тя Г (xr (1/2 1 і+ј+р 2 
Ls = = x 4 1/2 | i3 Α jap m 4λ. Ig + 4À L, 


(p9) 


where 
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A T (i+j+pt1/2 
+ј+р = п 


with 
Ау = 1 
amd 
Ал ΡΕΊ m ырыга Ass jap 
_ (itm a (x ,a)(-p,q) τ ш | x: °1/T 
2) (-1 " =: ο (luna) (тт) 79 p(1/2,x;x+1/2+a;1/x) 
(D10) 
Let: | 


Sp 
x x-m m l/2-m P(x)r(1/2 > 
R. = 1 + 2 T c ΤΉ fp Bp T 
| p=0 


| | 
= 1 + Gm) „1/2-ш T(x) (1/2 J у 
R, = 1 + — Tv T0172 


J=L p=0 


Cr] 
> 
H3 
+ 
hj 


i m -ma / I (х+1/2) ΤΟ fitjtp P 


i=l j=l p=O iJ 
(D11) 
Then equations (B38) for C and D can be written as follows: 
"TN Aik + 
C(A,m,k) = - ik + επ.) 
| 1 
А (D12) 
y? 2ik - 2k 1 Rz 
D = ik C - — - 


2 Ry 
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.076178 
. 065934 
050747 
. 040292 
.032811 
. 027299 
. 025138 
‚012545 


МАСА TN 5265 


TS82 


NACA TN 3263 


"= 
0 1 
1 1 
2 1 
5 1 
4 1 
S 1 

1 
1 
1 
1 
1 


TABLE II. 


39 


e 
Cute 


оошо мыо 


= 


a mæ 一 -~ 


0.008232 
. 007546 
. 006950 
«006429 
. 005970 
. 005563 
. 005200 
.004875 
.004582 
. 004318 
. 004078 
. 005860 
. 003660 
. 0053477 
. 003509 


а tg pe te a - 一 —  A— — . хш — — ‘M 


20 NACA TN 3263 


y / | 

— A (x, Yu 

/ W; Жы 
L y „7 | 
К dv. / 
/ ΕΓ 
РА / 
/ / 
/ 
/ Р, 
/ / 
/ / 


Figure 1. - Cascade geometry. 
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Figure 2. - Variation of Pj and Po with spacing. 
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Figure 3. - Variation of real and imaginary parts of C(A,k) with reduced 
frequency k for various values of spacing s. 
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Figure 4. - Variation of bending 
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Figure 5. - Variation of bending coefficient 
Reduced frequency, 0.4; phase lag m, 
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Figure 6. - Airfoil displacements. 
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Figure 7. - Contours in C-plane. 
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